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Abstract. The Yang-Lee, Fisher and Potts zeros of the one-dimensional Q-state 
Potts model are studied using the theory of dynamical systems. An exact recurrence 
relation for the partition function is derived. It is shown that zeros of the partition 
function may be associated with neutral fixed points of the recurrence relation. Further, 
a general equation for zeros of the partition function is found and a classification of 
the Yang-Lee, Fisher and Potts zeros is given. It is shown that the Fisher zeros 
in a nonzero magnetic field are located on several lines in the complex temperature 
plane and that the number of these lines depends on the value of the magnetic field. 
Analytical expressions for the densities of the Yang-Lee, Fisher and Potts zeros are 
derived. It is shown that densities of all types of zeros of the partition function are 
singular at the edge singularity points with the same critical exponent a = — i . 
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1. Introduction 

It is well established that the thermodynamic properties of a physical system can be 
derived from a knowledge of the partition function. Since the discovery of statistical 
mechanics it has been a central theme to understand the mechanism how the analytic 
partition function for a finite-size system acquires a singularity in the thermodynamic 
limit if the system undergoes a phase transition. The answer to this quest was given 
in 1952 by Lee and Yang in their famous papers [0. They considered the partition 
function of the Ising model as a polynomial in activity (exp(—2H/kT), where H is a 
magnetic field) and studied the distribution of zeros of the partition function in the 
complex activity plane. It was shown that phase transitions occur in the systems 
where continuous distribution of zeros of the partition function cuts the real axis in 
the thermodynamic limit. Also proved, was the circle theorem, which states that zeros 
of the partition function of the ferromagnetic Ising model lie on the unit circle in the 
complex activity plane ( Yang-Lee zeros). Later, Fisher initiated a study of zeros of the 
partition function in the complex temperature plane (Fisher zeros). Fisher showed that 
complex temperature zeros of the partition function of the Ising model in zero magnetic 
field on a square lattice lie on two circles \v ± 1| = ^/2, where v = tanh( J/2/cT). Since 
that time there has been a considerable body of work studying the partition function 
zeros of the Ising and Potts models on various regular lattices, most notably in recent 
years [|J. Zeros of the partition function were also studied for spin models defined on 
hierarchical || and recursive [|I(J lattices, random graphs [^TJ and aperiodic systems [12|] . 



spin glasses [Oj, percolation and Self-Organized Criticality models 

Nowadays, the investigation of zeros of the partition function becomes a powerful 
tool for studying phase transitions and critical phenomena. Recently, much attention has 
been attracted to the study of zeros of the partition function in one-dimensional systems 



in connection with helix-coil transitions in biological macromolecules |L5|. Hence, it is 
very important to study the general properties of zeros of the partition functions for 
different one-dimensional systems. 



In 1994, Glumac and Uzelac |T(| using the transfer matrix method studied the 



Yang-Lee zeros of the one-dimensional ferromagnetic Potts model for non-integer values 
of Q > 0. They showed that for < Q < 1 the Yang-Lee zeros are located on a 
real interval and for low temperatures these are located partially on the real axis and 
in complex conjugate pairs on the activity plane. Later on, Monroe [17) numerically 



studied the Yang-Lee zeros of the Potts model for some particular values of Q. Then, 



Kim and Creswick |T8[ showed that for Q > 1 the Yang-Lee zeros lie on a circle with 
radius R, where R < 1 for 1 < Q < 2, R > 1 for Q > 2 and R = 2 for Q = 2. 
Only recently the full picture of Yang-Lee zeros of the ferromagnetic Potts model for 



< Q < 1 has been found ]TD| . 

In this paper, the dynamical systems theory is used to study the Yang-Lee, Fisher 
and Potts zeros|| of the partition function for the one-dimensional Q-state Potts model. 

X Zeros of the partition function considered as a function of complex Q. 
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In Section 2, a recurrence relation for the partition function is derived. It is shown 
that zeros of the partition function may be associated with neutral fixed points of the 
recurrence relation. A general equation for zeros of the partition function is derived. 
Formulae for the free energy and the density of zeros of the partition function are found. 
In Sections 3 and 4 the method developed in Section 2 used to study the Yang-Lee and 
Potts zeros of ferromagnetic and antiferromagnetic Potts models. A classification of 
Yang-Lee and Potts zeros is given. In Section 5, the Fisher zeros in a nonzero magnetic 
field are investigated. It is shown that the Fisher zeros in a nonzero magnetic field 
located on several lines in the complex temperature plane and the density of Fisher 
zeros is singular at the edge singularity points with the same critical exponent as that 
of both the Yang-Lee and Potts zeros, a — — |. 

2. Zeros of the partition function of the Potts model 

The Hamiltonian of the one- dimensional Q-state Potts model in a magnetic field is 
defined as follows 



where 5 is the Kroneker delta function, <Tj denotes the Potts variable at site i and takes 
the values 0, 1, 2, . . . , Q — 1. The first sum in the r.h.s. of ([!]) goes over all edges and 
the second one over all sites on the lattice. For J > the model is ferromagnetic and 
for J < is antiferromagnetic. Note that due to the symmetry, the Hamiltonian ([]]) is 
the same if the external field H is applied to any spin state, namely, if 5(o~i, 0) in (0) is 
replaced by <5(<7j, a) for any a — 1, 2, . . . , Q — 1. For Q = 2 the Potts model corresponds 
to the Ising model and in order to keep the analogy with the Ising model we designate 
H as a magnetic field. We may assume the cyclic boundary condition a n = cr_ n and 
that the number of sites is 2n + 1 without losing a generality. Cutting the lattice at the 
central site do will separate it into two branches I and II with equal statistical weights 
g n (cr ) (Figure [I]). 



(i) 



<ij> i 




Figure 1. The procedure of derivation of the recurrence relation for the partition 
function. 
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Cutting the branch I (II) at the site o\ (c-i), the recurrence relation for g n (o~) may 
be found as 

Q-i 

9n(<ro) = ex P [^(°"o, en) + hS(ai, 0)] n _i(0i). (2) 

cti=0 

where J = J/kT and h = H/kT. The partition function may be written in the form 

z = J2e- H/kT = Eexp[WK0)]^( ffo ). (3) 

<T O"o=0 

Introducing the notation 

n <7n(<r = 0)' lj 
and using (0), the recurrence relation for x n may be found as 

Xn = /(*n-l), f(x) = Zfi+{Q _ 1)X , (5) 

where \i = e h , and does not have a direct physical meaning, but the 

thermodynamic functions, such as the magnetization, the specific heat and etc., may 
be expressed in terms of x n . For example, the magnetization for our model in the 
thermodynamic limit has the form 

m =( W Z)-.E,( ff ,0).-W = _^_, ( 6 ) 

where x is an attracting fixed point of the mapping (^|). It will be shown below how 
the thermodynamic properties of the model may be defined from the dynamics of the 
recurrence relation (H). 

The mapping (|5]) is a Mobius transformation, i.e. a rational map of the form 

R(x) = ^— , ad — be 0, 

cx + a 

where 

R(oo) = a/c, R(—d/c) = oo, 

if c 7^ 0, while R(oo) = oo when c = 0. The dynamics of such maps is well studied |H| . 
The mapping @ has only two fixed points which are solutions to the equation f(x) = x. 
According to the theory of complex dynamical systems, fixed points are classified as 
follows: A fixed point x* is attracting if |A| < 1, repelling if |A| > 1, and neutral if 
|A| = 1, where A = j-f(x)\ x=x * = f'( x *) is called the eigenvalue of x*. It is easy to 
show that either both fixed points of the mapping (El) are neutral, or that one of them is 
attracting while the other is repelling. The correspondence between thermo dynamical 
properties of the model and the dynamics of the mapping (|5]) is the following: If for a 
given temperature and magnetic field (z and fx) the mapping (|5|) has an attracting fixed 
point, then the system is in a stable state and its thermodynamical functions are defined 
by this fixed point (see for example @). The other fixed point is repelling and does not 
correspond to any phase. On the other hand, if the iterations of (EJ) do not converge to 
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a fixed point, i.e. the mapping (|5|) has neutral fixed points only, the system undergoes 
a phase transition. Therefore, the existence of neutral fixed points of (||) corresponds 
to a phase transition in the model. It was mentioned in the Introduction that zeros 
of the partition function correspond to phase transitions in the model. Thus, zeros 
of the partition function are associated with neutral fixed points of the corresponding 
mapping. These may be found from the conditions of existence of neutral fixed points 
of the mapping (|5|). These conditions are the following 

I = X (7) 

Excluding x from the system (|7|) after some algebra the equation of phase transitions 
may be found 

*V~2 [( z _ 1) ( z + Q _ 1) COS + 1_Q] fi + ( 2 + Q_ 2 ) 2 = 0, (8) 

where G [0, 2tt]. Solutions to this equation for different values of correspond to zeros 
of the partition function. Hence, the free energy of the model may be written in the 
form 

F ~ f 27T \n(z 2 fx 2 + 2 (Q - 1) fx + (z + Q - 2) 2 - 2 \i (z - 1) (z + Q - 1) cos ) d<f>. (9) 
Jo 

Formula @ may also be derived from the transfer matrix method [T^, ^TJ using two 



non-degenerate eigenvalues of the transfer matrix |2T 
1 



K,2 = \ zfi + z + Q-2±^/(zfi-z-Q + 2) 2 + A(Q-l)/2 
and the following mathematical identity for any pair of scalars C, D 
C N + D N = [] [C + exp(27rm/A^) D] , 



(10) 



where the product is from n = 1, 2, . . . , iV if N is odd; and from n = 1/2, 3/2, . . . , N — 
1/2 for N even. Moreover, one can easily show that the eigenvalues of fixed points of 
(|) are related to the eigenvalues of transfer matrix ([H]) as follows 

4fji{z-l)(z + Q-l) 



Ai, 



n,2) 

Hence, the condition of a phase transition |Ai| = [ A2 j (both Ai and A2 are neutral 
fixed points) corresponds exactly to the |A*| = jA^I condition for the eigenvalues of 
transfer matrix fliCf ). It means that the phase transition point based on dynamical 
systems approach coincides exactly with the phase transition point based on free energy 
considerations. This correspondence seems to be general for spin models defined on 
recursive lattices 0, |22|] . 

The density of zeros of the partition function may be found from the equation (|8|) 
in a standard way and has the form 

BdtA-AdtB 
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where = J^, 

A = z 2 fj 2 + 2 (Q - 1) fj. + (z + Q - 2) 2 , B = 2 /i (z - 1) (z + Q - 1) cos 0, 

and £ = fi, z or Q depending on whether the Yang-Lee, Fisher or Potts zeros are 
considered. It is interesting to note that in our case the numerator of (|TT1) always 
contains the multiplier (A+B). Hence, the density <?(£) is singular only when A—B = 0||. 
From (|TID follows that g(£) has a singular behavior g(£) ~ |£ — £*| CT , where is a solution 
to the equation A — B = 0. We will see that o~ = — \ for all types of zeros of the partition 
function. Also, note that the equation A — B = corresponds to the phase transitions 
equation (j8|) for <fi — 0. It defines the edge singularity points j2~0fl . In the subsequent 
Sections we will apply equations (||) and (0) to the study of the Yang-Lee, Potts and 
Fisher zeros. 



3. The Yang-Lee zeros 



According to the results of the previous section the Yang-Lee zeros of the Q-state Potts 
model may be found by solving the equation (|8|) with respect to \i. Equation @ is 
a quadratic equation of /i with real coefficients. Note that solutions to the equation 
(|) lie either on the real axis or in complex conjugate pairs on a circle with radius 
R = \z + Q — 2\/z and have the form 



where 



^1,2 



E 



E 



2 cos 2 f - F ± 2 



COS' 



I (cos- 



(z-l)(z + Q-V 



and F 



£ — p 



z{z + Q-2) 



(12) 



(z-l)(z + Q-l) 
with respect to \i has already been performed 



[TCfl . Here, we 



A detailed study of 

will present only the main results for the Yang-Lee zeros of both ferromagnetic and 
antiferromagnetic Potts models. 

For the ferromagnetic Potts model [z > 1) one can find that for Q > 1 all solutions 
([[ID are complex conjugate and lie on an arc of circle with radius R — (z + Q — 2)/z. 
Writing \i in the exponential form /i = Re %e the angular distribution of the Yang-Lee 
zeros may be found in the form 



e 

cos - 
2 



\ 



[z-i)(z + Q-r 



cos — . 
2 



(13) 



z(z + Q-2) 

From (|13|) one can see a gap in the distribution of Yang-Lee zeros, i.e. there are no 
solutions to the equation @ in the interval — 9 < 9 < 9 , where 9q = 2 arccos 
This is the well known gap in the distribution of Yang-Lee zeros of ferromagnetic models 
above the critical temperature (formally T c = for the one-dimensional case) first 
studied by Kortman, Griffiths and Fisher ||20||. The end points of the gap are called 



§ Solutions to the equation B=0 should be neglected since these do not correspond to zeros of the 
partition function. 
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Yang-Lee edge singularity points. From (^) and fll3|) follows that the Yang-Lee edge 
singularity points correspond to <fi = 0. Hence, according to ( |TTD the density of Yang-Lee 
zeros is singular in the Yang-Lee edge singularity points. Substituting = into ([E|) 
the formula for Yang-Lee edge singularity points we have 

fi± = ^[y/(z-l)(z + Q-l)± /T^g} 2 . (14) 

Substituting £ = /x = i?e l6i in (|TTJ) after some algebra, the density of Yang-Lee zeros for 
Q > 1 may be found in the form 

= l_ | sin f | 

27r */sin 2 § -sin 2 ^ 



where 9 = 2 arccos V From the equation (|I5|) follows that the density g(9) diverges 
in the Yang-Lee edge singularity points /i± with the critical exponent a = — ~, i.e. 
(7(6') oc \9 — 9 \~2 when <\> — > or — > O - 

For Q < 1 the Yang-Lee edge singularity points /i± are real and the density of 
Yang-Lee zeros has the form 



M = — ; '^-v^U (16) 
^^yl(ji + -y){ t i-H.) 

g(fi) diverges in the points /i±, i.e. g(/i) oc — /i±|°", with the critical exponent a = — |. 
The summary of results for the Yang-Lee zeros of the ferromagnetic Potts model is given 
in Figure 0. 




Figure 2. A schematic presentation of the Yang-Lee zeros of ID ferromagnetic Potts 
model. Here R = z+ ®~ 2 and fi± are defined in (|l4| ) . 



The Yang-Lee zeros of the antiferromagnetic Potts model (z < 1) may be studied 
in the same manner and the results are given in Figure |3|. For complex /i± the angular 
distribution of Yang-Lee zeros has the form 



9 

sin- 



(z-l)(z + Q- 1) 
\ z{z + Q-2) 



cos-. 



(17) 
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In contrast to the ferromagnetic case now the Yang-Lee zeros lie in the interval 
—9 <8<6 , where 9$ = 2 arcsin \/F~ l . The density of Yang-Lee zeros in this case has 
the form 

1 



9(0) 



cos • 



27r - /sin 2 2a - sin 2 2 



;is) 



For real values of /i± the density of Potts zeros has the same form as for the ferromagnetic 
case (1161) . Actually, the formula (III 



is the most general from which formulae (0) and 
([TJD may be derived (See also Section 4). Note, that here also the density of Yang-Lee 
zeros diverges at the Yang-Lee edge singularity points with the same index a = — |. 
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Figure 3. A schematic presentation of the Yang-Lee zeros of ID antiferromagnetic 
Potts model. Here R = 2 ~ ~ z and fi± are defined in (|i4|). 



4. The Potts zeros 

Recently, much attention has been attracted to the study of zeros of the partition 
function in the complex Q plane for the Potts model on regular lattices and finite 
graphs |[23|| . The partition function of the Potts model (H = 0) on a finite graph G, 
Zg{Q,v), may be considered as a polynomial in the number Q of Potts states and 
the temperature-like variable v = exp(J/fcT) — 1. At zero-temperature the partition 
function of the antiferromagnetic Potts model (v = —1) correspond to the chromatic 
polynomial Pg{Q), which is closely related to the Q-coloring problem. By definition, 
for graph G and positive Q, Pg{Q) is the number of ways in which the vertices of G can 
be assigned "colors" from the set 1, 2, . . . , Q in such a way that adjacent vertices always 
receive different colors. The original hope was that study of the real or complex zeros 
of Pg{Q) might lead to an analytic proof of the Four-Color Conjecture, which states 
that -Pg(4) > for all loop-less planar graphs. To date this hope has not been realized, 
although combinatoric proofs of the Four-Color Theorem have been found [Q. Even 
so, the Potts zeros and zeros of Pg{Q) are interesting in their own right and have been 
extensively studied in recent years . 

In this section the Potts zeros of ferromagnetic and antiferromagnetic one- 
dimensional Potts models are studied. It was shown in Section 2 that zeros of the 
partition function correspond to solutions of the equation (H), which is a polynomial 
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in Q. Hence, the Potts zeros may be found as solutions to (0) with respect to the 
parameter Q. It is convenient to rewrite (§) using the variable P = Q + z — 1 

P 2 -2P[l-/i + /i(z-l) cos0] + (l-^/i) 2 = O. (19) 

Equation ([OJ) is a quadratic equation in P with real coefficients. Solutions of ( ]T9| ) lie 
either on the real axis or on the circle with radius R = |1 — z \x\ in the complex plane 
P, and have the form 

P 1)2 = 1 - // + //(* - 1) cos0±2 cos | \ffi(z- 1)(1 -n-n(z- l)sin 2 |). (20) 

Since the analysis of ( ^0|) is standard, we will skip the details and give the results 
in Figures |] and |5] for ferromagnetic and antiferromagnetic Potts models respectively. 
In analogy with the Yang-Lee zeros, the Potts edge singularity points are defined as 
solutions to the equation (0) for = 0, and have the form 

Q± = P± ~ z + 1 = 1 - z + (y/l^i ± 2 (21) 
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Figure 4. A schematic presentation of the Potts zeros of ID ferromagnetic Potts 
model. Here R = z /J, — 1 and Q± are defined in (pl|). 
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Figure 5. A schematic presentation of the Potts zeros of ID antiferromagnetic Potts 
model. Here R = 1 — z fi and Q± are defined in 

The general formula of the density of Potts zeros has the from 

9(P) = ^-5 7=^=^ =,- ( 22 ) 
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From (p^ ) it follows that the density of Potts zeros is singular at the Potts edge 
singularity points. For complex values of Q±, (^) may be written in the form 

g[ff) = — ; |Si11 ^ = , (23) 

27r v / sin 2|_ sin 2|' 

for the ferromagnetic model, where \6\ > 9 and 9 = 2 arccos(/i (z — l)/(fiz — 
For the antiferromagentic model, g(6) has the form 

2 t ^/sin 2 f - sin 2 | 

where \Q\ < 9 and 9$ = 2 arcsin(/x(l — z)/(l — fiz))~2. Different formulae for the 
density of Potts zeros occur because the square root function in (E^) is not a unique 
function in the complex plane P, i.e. one branch of it corresponds to the ferromagnetic 
model and the other to that of the antiferromagnetic. 

It is noteworthy that the density of Potts zeros diverges at the Potts edge singularity 
points with the same critical exponent a = — | as for the Yang-Lee zeros. 

5. The Fisher zeros 

In this section the Fisher zeros in a nonzero magnetic field are studied. Usually, 
the Fisher zeros are considered as zeros of the partition function with respect to a 
temperature-dependent parameter. In our case z is such a parameter. Formally, the 
magnetic field in the Hamiltonian (HD may be presented in the form H = H J, where H 
is the renormalized magnetic field. Later we will refer to if as a magnetic field. Then, 
\x may be written as z H and the equation ([8]) will has the form 

P^(z,H,Q)=0, (25) 

where 

P^z, H, Q) = z 2H+2 - 2 cos0 z H+2 - 2(Q - 2) cos z H+l + 2(Q - l)z H + (z + Q - 2) 2 . 

P<j,(z, H,Q) is obviously a polynomial for integer values of H and (|2~5D may be solved 
numerically. For non-integer values of H the equation ( P5|) becomes a transcendental 
equation and the Fisher zeros may be found by numerically checking the condition of 
existence of neutral fixed points of the mapping @. 

The density of Fisher zeros may be found from (|TT|) by substituting £ = z and has 
the form 

( s 1 HG 1 (z,H,Q)-zG 2 (z,H,Q) 
g(z) = , , (2b) 

(z-l)z(z + Q-l) y/-P (z, H, Q) 

where 

Gi(z, H,Q) = {z-l){z + Q- l){z H+1 -z-Q + 2), 
G 2 (z,H,Q) = (z H - l)(2(z + Q - I) - Q z) + Q 2 . 
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From fl2"6|) it follows that the density of Fisher zeros is singular in the Fisher edge 
singularity points, which are defined as solutions of the equation ( p5|) for = 0. 

Let us first study the Fisher zeros for the Ising model (Q = 2). For Q = 2 the 
equation (|25| ) has the form 

z 2H+2 - 2 cos0 z H+2 + 4 cos 2 f z H + z 2 = 0. (27) 

This equation is symmetric under the H — > —H transformation because of the Z(2) 
symmetry of the Ising model. Since all coefficients of ([27|) are real, its solutions will be 
either real or complex conjugate. For integer values of H the Fisher zeros have the form 
shown in Figure |6|. From Figure |6] and (|27j) one can see that the Fisher zeros located on 




Figure 6. The Fisher zeros of ID Ising model for different values of a magnetic field. 
The big dots at the ends of the lines show Fisher edge singularity points. 



lines ending at the Fisher edge singularity points and the number of lines equals to the 
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value of the magnetic field \H\ + 1. Fisher edge singularity points are divided on pairs 
and every pair corresponds to a line. Note that z = is a twice degenerate Fisher edge 
singularity point with the critical exponent a = —\. Formally, two degenerate z = 
edge singularity points are considered here as forming a "line" . One can see that for odd 
values of \H\ there are negative Fisher zeros and for \H\ = 4n, n £ N, the Fisher zeros 
are partially located on the imaginary axis. Since P$(z, H, Q) is an analytic function of 
H the location of Fisher zeros for non-integer values of H are continuous deformations 
of Fisher zeros in the field [\H\] or [\H\] + 1, where [\H\] is the integer part of \H\, 
i.e. the minimal integer less than \H\. The location of Fisher zeros for non-integer H 
may be found numerically from the condition of existence of neutral fixed points of the 
mapping (|5|). As an illustration the Fisher zeros for non- integer values of < H < 2 
are given in Figure |f[ 

The Fisher zeros for Q ^ 2 Potts model may be studied in the same way as for the 
Ising model. It is impossible to give all possible configurations of Fisher zeros for any 
H and Q. Here we give only the summary of the main properties. First of all, there is 
no H — > —H symmetry for Q ^ 2 and the Fisher zeros are different for positive and 
negative values of H. 

For Q = 1 the Fisher zeros located on a closed curve and the density of Fisher zeros 
is not singular on this curve. 

For Q 1 Fisher zeros are located on several lines and the density of Fisher zeros 
is singular at the end points of these lines with the edge singularity exponent a — — |. 

Numerical experiments show that for integer values of H there are the following 
properties: for < Q < 1 there is an interval of real Fisher zeros only for even values 
of H; for 1 < Q < 2, the Fisher zeros intersect the negative real semi-axis only for odd 
values of H; for Q > 2 there is an interval of real Fisher zeros only for odd values of H. 
The number of the lines of Fisher zeros are defined as in the Ising model. In Figure [| 
we give some plots that illustrate these properties. 

6. Conclusions 

In this paper, the Yang-Lee, Fisher and Potts zeros of the one-dimensional Q-state Potts 
model are studied using the dynamical systems theory. A recurrence relation for the 
partition function is derived. It is shown that zeros of the partition function may be 
associated with neutral fixed points of the corresponding recurrence relation. A general 
equation for zeros of the partition function is derived. It is shown that the density of 
zeros of the partition function is singular in the edge singularity points with the critical 
exponent o = — |. 

The recursive method used in this paper is equivalent to the transfer matrix method, 
but in some sense it is more flexible. It gives a numerical algorithm to study the zeros 
of the partition function. The algorithm is based on testing the condition of existence 
of neutral fixed points of the corresponding recurrence relation: If neutral fixed points 
exist, it means that the partition function is zero for given values of z, fi and Q. In 
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Figure 7. The Fisher zeros of ID Ising model for non-integer values of a magnetic 
field 0<H<2. 



Section 5 we saw that this algorithm is very useful for studying the Fisher zeros in an 
arbitrary magnetic field. Moreover, it is the only method for studying the Fisher zeros 
for non-integer values of if, where H is the renormalized magnetic field H and H — H/J. 
In conclusion we would like to note that the results given in this paper show that the 



Partition function zeros of the one- dimensional Potts model 



14 




thermodynamic properties of the one-dimensional Potts model are completely defined 
by the recurrence relation (^) or equivalently, by two non-degenerate eigenvalues of 
the corresponding transfer matrix (|10|) (see also fllQl). Moreover, it is proved that the 



phase transition point based on dynamical systems approach coincides exactly with 
the phase transition point based on free energy considerations. The dynamical systems 
approach used here gives a numerical method which may be used for studying other one- 
dimensional systems for which a one- dimensional recurrence relation may be derived. 
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